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VAPORIZATION  WAVE  TRANSITIONS* 


ABSTRACT 

A  vaporization  wave  analogous  to  the  unsteady  simple  wave  in  an 
expanding  gas  is  proposed  to  account  for  phenomena  observed  in  the 
evaporation  of  superheated  metals.  The  vaporizing  model  is  visualized 
on  thermodynamic  grounds  as  carrying  the  liquid  metal  through  a 
continuous  succession  of  states  either  on  or  near  the  liquidus  line  in 
the  two-phase  region.  On  this  line,  the  adiabatic  sound  speed  for  wet 
vapor  will  limit  the  rate  of  propagation  of  the  vaporization  front 


*This  report  contains  the  essence  of  two  lectures  given  oy  tne  autnor 
to  the  International  School  of  Physics  " Enrico  Fermi"  Course  on 
Physics  of  High  Energy  Density t  1969 ,  Varenna,  Italy. 
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into  the  liquid.  Experimental  data  for  wire  explosions  of  Al,  Ag,  Cu, 
Au,  Fb  and  Hg  (frozen)  are  analyzed  for  wave  speeds.  While  the  in¬ 
fluence  of  thermal  expansion  of  the  liquid  can  be  accounted  for  the¬ 
oretically,  insufficient  thermal  data  are  available  for  the  metals  to 
permit  correction  of  the  wave  speeds  for  this  effect.  The  experimen¬ 
tally  derived  wave  speeds  are  compared  with  theoretical  values  of  the 
adiabatic  sound  speed  in  the  wet  vapor  obtained  from  a  modified  van 
der  Waals  equation  of  state.  At  low  velocities,  the  agreement  is 
satisfactory  but  higher  values  deviate  considerably  from  theory. 
Possible  causes  of  the  deviations  include  the  crudity  of  the  fluid 
dynamic  model,  neglect  of  thermal  expansion,  lack  of  information  about 
the  relationship  between  density  and  electrical  conductivity  and  the 
approximation  imposed  by  the  van  der  Waals  equation. 
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1.  INTRODUCTION 


We  outline  here  some  theoretical  and  experimental  aspects  of  the 
dynamical  process  by  which  an  impulsively  heated  fluid  makes  the 
transition  from  liquid  to  vapor.  Our  interest  is  in  the  initial 
stages  of  the  expansion  as  it  occurs  in  the  vaporization  of  super¬ 
heated  metals. 

For  polytropic,  non-reacting  fluid  the  head  of  an  expansion  wave 

proceeds  into  the  compressed  region  at  the  local  velocity  of  sound  and 

1* 

the  ensuing  expansion  is  a  simple  wave.  When  the  fluid  can  be  con¬ 
densed,  the  first  expansion  will  also  be  limited  by  the  speed  of  sound 
in  the  liquid;  but  the  vaporization  which  follows  may  take  place 
either  from  the  bounding,  external  surface  or  by  cavitation  with  vapor 
bubble  formation.  Since  formation  of  internal  cavities  is  limited  by 
inertia  and  nonequilibrium  effects,  and  may  require  a  considerable 
time  interval  for  bubble  formation,  we  shall  consider  only  experimental 
situations  where  vaporization  takes  place  from  the  exposed  surface. 

In  this  case  the  evaporation  will  be  limited  by  some  maximum  speed  with 
which  the  head  of  the  wave  travels.  Because  a  phase  change  occurs  and 
because  the  resulting  damp  vapor  is  more  compressible  than  the  liquid, 
one  would  expect  the  limiting  wave  speed  to  be  much  slower  than  that  of 
sound  in  the  homogeneous  liquid.  In  other  words,  the  rate  at  which  a 
uniform  vaporization  phenomenon  should  occur  will  be  much  slower  than 
the  rate  at  which  a  purely  liquid  expansion  is  propagated. 


*References  are  listed  on  page  29. 
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2.  WAVE  HYPOTHESIS 


We  wish  to  test  the  hypothesis  that  the  limiting  upper  velocity 
of  a  dynamic  vaporization  process  will  be  a  sound  speed  characteristic 
of  the  two-phase  region  of  the  fluid.  In  particular  the  chosen  speed 
will  be  that  of  the  coexistence  side  of  the  liquidus  line  which  is  the 
low  volume  boundary  of  the  two-phase  region,  represented  as  the  locus 
Va(T)  in  Figure  1. 

For  a  substance  characterized  by  an  equation  of  state  P  =  P(v,  T), 
where  V  =  l/p  is  specific  volume  and  P,  p,  T  refer  to  pressure, 
density  and  temperature  respectively,  one  can  readily  show  that  the 
speed,  c,  of  a  small  disturbance  is  given  by  c2  =  (dP/dp)s  ,  or  equiv¬ 
alently  by 

c2  =  (SP/Sp)T  +  (t/ p2Cv  )  (dP/ST ) ^  ,  (l) 

where  Cy  is  the  specific  heat  at  constant  volume.  We  note  that  within 
the  two-phase  region  the  number  of  free  variables  is  reduced  to  one 
and  P  =  p(t),  the  vapor  pressure;  thus,  only  the  second  term  of 
Equation  (l)  remains  and, 

c  =  (T/Cvp2)^  (dP/dT)  .  (2) 

On  the  boundary  given  by  the  liquidus  line  two  values  of  sound  speed 
occur,  viz  ,  those  for  the  liquid  and  liquid  plus  damp  vapor.  Sound 
speed  is  double  valued  on  the  loci  where  slope  discontinuities  occur 
in  the  adiabats . 

We  refine  the  wave  hypothesis  slightly,  to  state:  dynamic 
vaporization  phenomena  in  superheated  liquids  move  into  the  undisturbed 
fluid  no  more  rapidly  than  the  liquid-vapor  sound  speed  characteristic 
of  the  liquidus  line.  A  more  detailed  development  of  the  thermo¬ 
dynamical  model  will  be  given  in  §  6  below. 


10 


While  our  main  experience  with  the  vaporization  waves  comes  from 
experiments  with  metals,  one  would  expect  the  hypothesis  to  apply 
equally  well  to  a  broader  class  of  liquids  including  elements,  com¬ 
pounds,  solutions  and  mixtures.  Necessary  conditions  for  any  liquid 
to  exhibit  vaporization  wave  phenomena  can  be  tentatively  stated  as 
follows .  Such  a  substance  should  have  l)  a  piecewise  continuous 
equation  of  state  relating  pressure,  specific  volume  and  temperature, 

2)  a  change  in  phase  with  a  large  increase  in  specific  volume,  3)  a 
latent  heat  of  vaporization  large  compared  with  the  specific  energy 
at  the  boiling  point;  and,  for  practical  reasons,  4)  a  temperature 
high  enough  so  that  specific  energy  exceeds  a  certain  threshhoid. 

Wave-like  vaporization  phenomena  are  encountered  in  the  steam 
shock  tube.  Here,  water  contained  at  temperatures  and  pressures 
higher  than  its  normal  boiling  point  is  suddenly  released,  by  breaking 

a  diaphragm,  into  a  driven  section  containing  air  or  other  gas  at  lower 

"  ;  4 

pressure.  Experiments  or  this  kind  have  been  performed  by  Brown  and 

5 

by  Terner.  Brown  treats  the  expansion  wave  traveling  into  the  driver 
by  assuming  a  two-stage  process.  In  the  first  stage  the  water  over¬ 
expands  to  a  lower-pressure,  metastable  state.  In  the  second  it 
undergoes  an  exponential  relaxation  with  increase  of  specific  volume 
to  its  final  equilibrium  state.  The  connection  between  metastable  and 
equilibrium  states  is  found  by  using  certain  additional  assumptions 
and  the  conservation  laws  of  fluid  mechanics  in  their  njump!f  form. 

While  there  is  some  arbitrariness  about  choosing  the  metastable  state, 
Brown1 s  experimental  data  appear  to  agree  somewhat  better  with  the 
results  of  the  two-stage  theory  than  with  those  of  an  equivalent 
equilibrium  expansion. 
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3-  EXPLODING  WIRE  EXPERIMENTS 

From  the  point  of  view  of  historical  order,  the  vaporization  wave 
hypothesis  was  first  developed  to  account  for  some  of  the  peculiar 

phenomena  observed  when  wires  are  exploded  by  a  condenser  discharge. 

6 

Somewhat  later  it  was  proposed  as  a  more  general  concept  applying  to 
the  explosive  expansion  of  superheated  liquids  in  general. 

In  an  exploding  wire  experiment  a  long,  thin,  conducting  cylinder 

is  impulsively  heated  by  an  energetic  condenser  discharge.  Depending 

on  the  choice  of  the  conductor  and  the  initial  conditions  a  variety  of 

complex  phenomena  occur  in  the  ensuing  expansion.  We  shall  not 

describe  the  entire  process  here  but  for  a  more  complete  discussion 

7 

refer  the  reader  to  a  recent  review.  What  will  interest  us  are  the 
early  stages  of  the  explosion  as  depicted  in  Figure  2.  There  one  sees 
a  streak  camera  record  of  the  expansion  of  a  cylinder  of  copper  cor¬ 
related  in  time  with  simultaneous  oscilloscope  traces  of  the  current 
passing  through  the  wire  and  the  voltage  measured  across  it.  The 
small  ramp  in  the  voltage  curve  represents  melting  of  the  wire. 

Shortly  thereafter  the  streak  photograph  shows  a  haze  of  vapor  expand¬ 
ing  about  the  more-dense,  interior  of  the  wire.  During  this  process 
the  current  decays  rapidly  to  zero  in  the  same  interval  that  voltage 
goes  through  a  peak,  which  in  this  case  is  about  four  times  as  high 
as  the  initial  condenser  voltage.  For  the  example  chosen,  the  wire  is 
so  well  matched  to  the  condenser  that  nearly  the  entire  initial  energy 
is  deposited  during  the  first  pulse.  The  physical  evidence  for  a 
vaporization  wave  is  contained  in  the  interval  of  wedge-shaped  expan¬ 
sion  terminated  by  the  zero  of  the  current  pulse. 

An  interpretation  of  the  exploding  wire  data  can  be  obtained  as 
follows.  To  first  approximation  we  suppose  the  wire  material  to 
exist  in  two  states  only,  viz.,  (l)  fully  conducting  metal  not  yet 
affected  by  the  expansion  wave,  and  (2)  expanded,  nonconducting,  wet 
vapor.  Thus,  for  a  vaporization  wave  propagating  with  velocity  v(t) 
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inward  from  the  wire  periphery ,  the  radius  of  the  conducting  core 

rt 

diminishes  with  time  as  r  =  r0  [l  -  (v(5)/r0)  dg ]  ,  where  t0  is 

-  to 

the  time  chosen  to  represent  the  start  of  the  expansion.  In  the 
example  of  Figure  2,  t0  would  be  a  little  more  than  4  p,sec.  With 
this  model,  conduction  will  cease  when  r  =  0.  Thus,  the  model 
provides  an  interpretation  of  the  current  shut-off;  and,  in  fact, 
initiates  an  explanation  of  dwell. 


4.  CONSTANT  VELOCITY  VAPORIZATION  WAVES 

Since  wire  resistance,  R,  varies  inversely  with  area  of  cross 
section  we  can  write 


(v(|)/r0)  d§] 


(3) 


If,  for  simplicity,  we  neglect  temperature  effects  and  assume  the  wave 
speed  to  be  constant ;  then,  with  t  =  r0/v  and  st  =  t  -  t0  , 

Equation  (3)  becomes 


3 

R  =  Ro/(l-s) 


(4) 


As  Figure  2  shows,  the  expansion  takes  place  in  a  short  time 
interval  near  the  maximum  of  current.  When  current  is  at  its  maximum 
the  charge  on  the  condenser  is  zero  and  the  energy  supply  is  stored 
in  the  magnetic  field.  If  we  assume  condenser  voltage  to  be  negligible 
during  the  current  decay,  the  circuit  equation  becomes 
L(di/dt)  +  Ri  =  0;  and,  with  the  previous  work,  can  be  written 

2 

di/ds  +  ai/(l-s)  =  0  ,  (5) 

where  a  =  (r0/v)(R0/L)  is  the  ratio  of  the  time  constant  of  the  wave 
to  that  of  the  circuit.  Integrating  Equation  (5)  yields 

i  =  i0  exp[ -  as/ (l-s) ]  ,  (6) 
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and  with  the  definitions  for  voltage,  VR  and  power,  PR  , 

Vr  =  R0 lo  (1_s)  exp[ -  as/(l-s)]  >  (7) 

and 

—2 

Pr  =  R0io  (l-s)  exp[-  2as/(l-s)].  .  (8) 

Differentiation  shows  that  VR  and  PR  have  maxima  at  1  -  a/2  and  1  -  a 
respectively.  Furthermore,  these  maxima  have  their  minimum  values,  as 
functions  of  a,  at  a  =  2  and  a  =  1  respectively.  A  variety  of  dif- 

O 

ferent  pulse  shapes  can  be  interpreted  with  this  elementary  calcula¬ 
tion.  Broadly  speaking,  if  the  time  for  e-fold  decay  of  current  is 
small  compared  with  the  time  of  wave  passage,  i.e.,  a  >  2,  the  current 
and  voltage  decay  curves  have  exponential  shapes  without  peaks.  Con¬ 
versely,  if  a  <  1  and  the  wave  is  comparatively  fast,  then  the  current 
decay  is  quite  flat  at  first  and  steepens  rapidly  as  the  wave  approaches 
the  center.  Voltage  passes  through  a  peak  whose  height  increases  as  a 
gets  smaller.  The  present  treatment  can  be  generalized  to  any  portion 
of  the  current  cycle  providing  only  that  the  vaporization  wave  relax¬ 
ation  time  is  much  smaller  than  a  quarter  period  of  the  undamped 
circuit  so  that  the  condenser  voltage  can  reasonably  be  taken  constant. 
The  analysis  of  the  constant  speed  vaporization  wave  shows  clearly  the 
interrelation  between  circuit  and  wave  speed  parameters;  and  shows, 
furthermore,  the  possibility  of  a  unified  discussion  of  a  wide  variety 
of  pulse  shapes.  It  fails  to  provide  an  accurate  method  of  calculating 
actual  wave  speeds. 


5 .  WAVE  SPEEDS  FROM  EXPERIMENTS 

If  one  abandons  the  program,  initiated  with  the  discussion  of 
constant  speed  waves,  of  attempting  to  provide  a  wave -speed  theory  of 
current  and  voltage  pulses,  a  considerable  shift  in  point  of  view 
becomes  possible.  We  regard  the  i-VR  curves  as  the  raw  data,  the 
physical  evidence  of  the  vaporization  wave,  and  attempt  to  deduce  wave 
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speeds  from  them.  Several  additional  assumptions  must  be  made  in 

order  to  provide  a  rational  reduction  scheme.  We  outline  it  here;  the 

3  9 

details  may  be  found  elsewhere.  5 

If  we  neglect  all  work  and  heat  terms  except  Ohmic  heating,  con¬ 
servation  of  energy  allows  us  to  equate  rise  in  specific  energy  of  the 
current-carrying  core  with  the  electrical  energy  deposited*  From  the 
corrected  voltage  VR  and  current  i  curves  one  can  calculate 


resistance  R  =  VR/i,  power  PR  =  VR i  and  energy  E 


VR i  dt .  We  can  , 


wr  ite 


de  =  i2  R  dt/m  ,  (9) 

=  V?  dt/Rm  ,  (10) 

where  e ,  m  are  respectively  the  specific  energy,  and  the  mass  heated 

during  interval  dt.  The  quantity  Rm  is  an  invariant  under  the 

2 

vaporization  wave.  The  factor  (r/r0)  cancels  out  of  the  product; 
however,  temperature  variation  remains  and  must  be  accounted  for. 

To  this  end  we  further  assume  that  resistance  is  linear  with 
specific  energy  above  the  melting  point.  This  assumption  can  be  ex¬ 
pressed  by  writing 

R  =  Ro  [1  +  P  (e-e0)](r0/r)2  ,  (ll) 

where  3  is  a  constant  to  be  determined  from  the  experiment.  If 
specific  heat  s  is  constant  in  the  interval,  then  3  =  cy/s,  e  =  sT  and 
Equation  (ll)  contains  the  usual  statement  that  resistance  increases 
linearly  with  temperature.  Using  Equations  (10)  and  (ll)  and  inte¬ 
grating,  we  find 

(e-e0)  +  (|3/2)(e-e0)2  =  (l/m^)  f  vj  dt  .  ( 12 ) 

J  t 

-o 


15 


To  determine  time  t0  when  the  vaporization  wave  begins,  one  plots 
scaled  resistance,  or  resistivity,  against  apparent  specific  energy 
E/m0  .  Before  vaporization,  the  radius  is  taken  to  be  r  =  r0  and  the 
constant  g  in  Equation  (ll)  can  be  determined  from  the  linear  portion 
of  the  resistivity  data  above  melt.  After  vaporization  commences,  R 
rises  steeply  away  from  the  linear  law  and  this  law  is  extrapolated 
to  account  for  the  effects  of  temperature  variation.  The  point  at 
which  R  departs  from  the  line  determines  R0  and  e0  .  Since  R  is  al¬ 
ready  known  as  a  function  of  time  t0  is  then  determined.  True  (e,t)  ■ 
values  may  be  obtained  by  integrating  Vp  as  in  Equation  ( 12 ) ;  and 

(e,R)  values  by  comparison  with  the  (R,t)  data.  With  (e,R)  data 

2 

Equation  (ll)  may  be  solved  for  (r/r0 )  .  The  resulting  points  are 
fitted  to  an  interpolation  curve  and  differentiated  numerically  to 
obtain  wave  speed.  The  initial  radius  r0  is  taken  nominally  to  be 
that  of  the  cold  wire.  A  more  elaborate  reduction  method  has  been  de- 

Q 

vised  to  correct  for  thermal  expansion  of  the  liquid  cylinder;  but, 
because  thermodynamic  data  for  metals  at  high  temperatures  are  lacking, 
has  not  yet  been  used.  The  correction  if  it  could  be  made  would  tend 
to  increase  the  values  of  wave  speed  deduced  from  measurement  by  a 
variable  factor  which  increases  approximately  from  1  at  room  tempera¬ 
ture  to  1.7  at  the  critical  temperature. 

The  assumption  of  a  linear  dependence  of  resistivity  on  specific 
energy  is  crucial  to  the  reduction.  Without  it  no  account  can  be  taken 
of  the  effects  of  temperature  rise,  and  these  are  known  to  be  too  large 
to  neglect.  Certain  metals,  e.g.,  Fe,  Ni,  W  show  anomalous  decreases 
of  resistance  above  the  melting  point  in  an  exploding  wire  experiment. 
For  this  reason  no  wave  speeds  have  yet  been  obtained  for  these  metals 
or  for  others  like  them.  If,  as  suspected,  the  cause  lies  in  early 
voltage  breakdown  and  associated  current  conduction  in  paths  external 
to  the  wire  there  is  some  hope  that  a  dense  ambient  fluid  may  prevent 
the  charge  leakage  and  enable  wave  speed  measurements  to  be  made. 


Wave  speeds  for  Cu  have  been  deduced  from  electrical  measurements 

as  described  above  and  compared  with  theoretical  values  calculated 

from  Equation  (2),  Values  of  liquid  density,  specific  heat  and  vapor 

pressure  were  obtained  from  tabulated  values  as  described  in  Reference 

6.  Owing  to  the  scarcity  of  measured  values  of  these  quantities  at 

the  desired  specific  energies,  various  approximations,  e.g.,  constancy 

of  liquid  density  and  specific  heat,  had  to  be  made.  Nevertheless, 

3 

agreement  between  measured  and  calculated  wave  speeds  for  Cu  ,  in  the 
temperature  range  from  the  boiling  point  to  the  critical  temperature, 
is  quite  satisfactory  and  a  considerable  encouragement  to  further  work. 
Accordingly,  similar  measurements  were  made  on  metals  A1  and  Fb.  The 
experimental  values  may  be  seen  in  Figure  3?  the  theoretical  curve 
shown  there  comes  from  a  thermodynamical  model  and  assumed  equation  of 
state,  see  §  6. 

6 

When  the  tabulated  data  were  used  to  calculate  theoretical  values 
for  wave  speeds  in  Fb ,  poor  agreement  with  experiment  was  obtained  and 
the  absence  of  reliable  high -temperature  vapor  pressure  data  for  this 
element  was  thought  to  be  responsible.  Since,  generally  speaking,  Hg 
is  the  only  metal  for  which  measurements  exist  of  pressure,  density 
and  conductivity  up  to  the  critical  temperature,  there  is  not  much 
hope  that  expected  values  of  vaporization  wave  speed  can  be  calculated 
for  other  metals  with  sufficient  precision  to  provide  a  correct  trend 
of  the  function  for  comparison  with  experimental  values. 

6.  THERMODYNAMIC  MODEL 

Accordingly,  we  develop  here  a  simplified  fluid  model  to  repre¬ 
sent  the  transient  behavior  of  a  material  heated  rapidly  from  the 
solid  state  up  through  its  critical  temperature.  Our  motivation  is 
to  supply  a  simple  yet  reasonably  accurate  basis  for  describing  the 
state  and  dynamical  behavior  of  superheated  liquids.  In  order  to  make 
any  further  detailed  flow  calculations  an  equation  of  state  is  needed; 
accordingly,  the  goal  will  be  to  find  an  approximate  equation  of  state 
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which  is  capable  of  representing  what  is  known  about  liquid  metals  and 
capable  of  filling  in  the  large  gaps  where  no  data  are  available. 
Agreement  will  be  sought  with  vaporization  wave  speeds,  for  example, 
even  though  detailed  internal  checks  with  other  quantities  are  impos¬ 
sible*  Our  model  is  based  on  the  following  hypotheses:  l)  the  tem¬ 
perature  rise  and  transitions  take  place  under  conditions  of  local 
thermodynamic  equilibrium,  2)  a  form  of  van  der  Waals  equation  ap¬ 
plies,  3)  realistic  specific  heats  can  be  assigned  the  liquid  and 
vapor  states  and  U)  the  liquid  expansion  takes  place  along  the 
liquidus  line,  V3(t)  of  Figure  1. 

The  first  assumption  is  one  mainly  of  simplification  although  a 

3 

plausible  argument  can  be  given  that  vaporization  from  the  free 
surface  should  occur  nearly  at  equilibrium  conditions. 

Van  der  Waals  equation  is  representative  of  a  class  of  state 
equations  which  are  reasonably  accurate,  though  not  precise,  and  yet 
embody  the  main  features  of  most  condensable  gases,  viz.,  a  hard  core 
repulsion  and  a  weak,  long-range  attractive  force.  With  the  modifica¬ 
tions  to  be  discussed  below  van  der  Waals  equation  allows  different 
and  realistic  specific  heats  to  be  assigned  the  liquid  and  vapor 
phases,  a  matter  which  improves  later  agreement  with  experiment. 

Choice  of  the  liquidus  line  as  the  thermodynamic  path  is  somewhat 
arbitrary  although  it  appears  to  be  reasonable,  as  may  be  seen  from 
the  following  argument.  For  most  physical  examples,  the  isotherms  and 
adiabats  in  the  liquid  phase  are  very  steep  compared  to  those  of  either 
the  coexistence  or  vapor  regions,  except  very  near  the  critical  point. 
Slope  magnitudes  of  the  adiabats  are  proportional  to  sound  speed 
squared;  consequently  the  sound  speed  in  the  liquid  is  much  higher 
than  in  either  of  the  other  regions.  We  may  then  argue  that  heating 
moves  the  thermodynamic  state  along  the  saturated  liquid  line.  Any 
compressive  tendency  to  drive  the  expanding  system  into  the  all-liquid 
state  will  be  rapidly  counteracted  by  liquid  thermal  expansion,  which 
is  dominated  by  the  fastest  speed  of  sound.  The  liquid  expansion 
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lowers  the  pressure  to  that  of  the  iiquidus  boundary  where  any  further 
expansion  must  be  accompanied ,  under  the  assumed  equilibrium  condi¬ 
tions.  by  partial  vaporization,  which  is  governed  by  the  slower 
vaporizing  wave  speed.  For  interior  parts  of  the  system  where  in¬ 
ertial  confinement  prevails,  premature  expansion  via  incipient  cavi¬ 
tation  will  cause  local  pressure  rise  which  will  force  the  fluid 
element  back  toward  the  Iiquidus  line.  This  process  of  excursion 
about  the  Iiquidus  continues  during  the  electrical  heating  pulse 
until  the  fluid  particle  evaporates  from  the  surface  by  expansion 
through  the  vaporization  wave. 

In  what  follows  we  take  the  thermodynamic  system  of  interest  to 
be  a  uniform,  molten  cylinder  of  metal  to  which  heat  is  being  added 
and  from  whose  surface  partial  vaporization  occurs.  The  speed  of  the 
leading  edge  of  the  vaporizing  wave  is  assumed  to  proceed  with  the 
wave  speed  of  the  saturated  wet  vapor  at  absolute  temperature  T,  as 
given  by  Equation  (2)  evaluated  on  the  Iiquidus  line. 

For  unit  mass  of  fluid,  pressure  P  and  specific  volume  V  are  re¬ 
lated  by  the  van  der  Waals  equation  in  the  form 

Pw  =  RT/ (V-b)  -  a/vS  >  (13) 

where  the  constants  a,  b  are  given  in  terms  of  values  at  the  critical 
point  by 

a  =  3PcVc  ,  b  =  Vc/3  ,  RTC  =  (8/3)  PCVC  .  (l4) 

Inequalities  T  >  0,  V  s  b  are  understood  to  hold.  When  T  ^  Tc  the 
van  der  Waals  equation  of  state  describes  a  two-phase  region  where 
both  liquid  and  vapor  phases  coexist,  as  in  Figure  1.  In  the  co¬ 
existence  region,  vapor  pressure  PA  is  not  derived  directly  from 
Equation  (13)  which  has  the  form  of  a  cubic,  but  rather  from  a  gen¬ 
eralization  of  the  Maxwell,  equal  area  rule  which  replaces  the  two 
loops  by  a  horizontal  line  cutting  equal  areas  from  the  loops. 
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G.  D.  KahlJ"'J  has  proven  the  remarkable  result  that  the  equal  area 
condition  is  not  necessary,  but  defines  only  one  of  a  continuum  of 
possible  choices.  For  T  <  Tc  he  writes 

rvi  (T) 

Pw  dV  =  pa(t)  [vx(t)  -  v3(t)]  +  0(t)  ,  (15) 

v3(t) 

where 

T 

0(T)  =  I  [Cv (t , liq)  -  Cv (t ,vap) ] (l  -  T/r)dT  .  (l6) 

w  T 

Specific  heat  at  constant  volume  Cv  is  a  function  not  only  of  temper¬ 
ature  but  also  the  state  of  the  fluid,  whether  liquid  or  vapor.  As 
Kahl  shows,  the  features  of  equilibrium  thermodynamics  are  preserved, 
the  main  one  being  the  stationary  property  of  the  Gibbs  function  as 
the  system  passes  from  liquid  to  vapor  on  an  isotherm.  If  Cv  is  the 
same  for  both  liquid  and  vapor,  then  0(t)  =0  and  Equation  (15)  gives 
the  Maxwell  result.  If  the  two  specific  heats  are  not  the  same,  then 
0(T)  /  0,  the  areas  are  not  equal  and  some  further,  useful  flexibility 
is  available  in  making  a  realistic  choice  of  the  two  specific  heat 
functions. 

To  evaluate  Equation  (2)  for  the  vaporization  sound  speed  on  the 
liquidus  line,  PA  and  V3  =  l/p3  can  be  found  numerically  from  Equations 
(15)  and  (l6).  The  specific  heat  for  the  coexistence  state,  CVA  , 
must  also  be  found  formally  and  then  numerically.  To  do  this  one 
finds  the  internal  energy  EA  and  differentiates,  using  CVA  =  (dEA/dT)v. 
Ea  may  be  determined  from  the  general  relation 

(SE/dV)T  =  T(dP/dT)v  -  P  ,  (17) 

and  noting  from  Equation  (15)  that  M(t)  =  T(dPA/dT)v  -  PA  is  independ¬ 
ent  of  V.  Direct  integration  of  Equation  (17)  gives 

EA  =  Ew  +  M(T)[V  -  V3(T)]  ,  (18) 
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where  Ew  =  Ew (v3 ,  T)  is  the  internal  energy  on  the  liquidus  given  by 
the  van  der  Waals  function.  One  finds 

CVA  =  (SEW/ST)V  +  (dM/dT)(V-V3)  -  M  (dV3/dT)  ,  (19) 

and 

(aEw/5T)v  =  (SEw/sV3)t  (dV3/dT)  +  (SEw/dT)v  ;  (20 ) 

finally,  on  the  saturated  liquid  line 

CVA  =  [(dEw/5V3)T  -  M]  (dV3/dT)  +  (5EW/ST)V  .  (21) 

The  last  terra  is  Cv(T,liq)  and  may  be  taken  from  experiment  or  repre- 
sented  approximately  in  some  chosen  functional  form.  Van  der  Waals 
equation  gives  no  further  information.  Putting  Pw  in  Equation  ( 17 ) 
gives  BEw/SV3  =  a/v|  .  With  these  results,  Equations  (15),  (l6),  (2l) 
and  van  der  Waals  equation,  the  vaporization  wave  speed  can  be  numer¬ 
ically  evaluated  on  V3(t).  In  the  liquid  the  more  general  definition 
of  Equation  (l)  must  be  used. 

While  the  equation  of  state  is  expressed  in  terms  of  absolute 
temperature,  the  independent  variable  one  obtains  more  easily  from 
experiment  is  the  added  specific  heat  content  Aq.  One  can  make  the 
connection  between  q  and  T  by  integrating  dq  =  dE  +  pdV  along  V3(t) 
from  the  melting  temperature  TM  to  T  using  computed  values  of 
Ew(Va,T)  and  Pa(t).  For  any  T  ^  Tc  a  correspondence  is  thus  given 
between  Aq(T)  and  T.  One  then  writes  q(T)  =  q(TM)  +  Aq(T)  where 
q(TM)  =  qM  is  the  heat  content  of  the  liquid  at  melt,  assumed  known 
from  other  sources.  Clearly  the  scaled  variable  q  =  q/RTc  depends 
upon  the  specific  heats  of  the  liquid  and  vapor  phases  through 
Equations  (15)  and  (l6),  as  well  as  on  qM/RTc •  This  latter  quantity 
has  a  value  near  0.60  for  a  number  of  metals.  Both  scaled  wave  speed 

i 

c  =  Cw/(RTc)  and  scaled  heat  content  are  functionals  of  the  specific 
heats  for  liquid  and  vapor.  When  one  can  approximate  both  of  the 
specific  heat  functions  by  constants,  the  parametric  dependence  of 
the  wave  speed  curve  on  the  choices  of  specific  heats  can  readily  be 
shown  as  in  Reference  9* 
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7.  COMPARISON  OF  THE  THERMODYNAMIC  MODEL  WITH  EXPERIMENT 

Wave  speed  data  for  six  metals,  Au,  Ag,  Cu,  Al,  Fb,  Hg,  are 

displayed  in  Figure  3*  While  a  somewhat  similar  plot  appears  in 

Reference  9  where  cross-hatched  areas  are  used  to  represent  the  more 

abundant  data,  here  the  individual  data  points  are  included  for  a 

range  of  variables  which  includes  most  of  the  two-phase  field. 

Dimensionless  wave  speeds  c  are  plotted  against  dimensionless 

specific  heat  content  q  .  Because  measured,  thermal  expansion  data • 

for  liquid  metals  are  lacking,  no  correction  could  be  made  for  the 

effects  expansion  prior  to  the  passage  of  the  vaporization  wave. 

Likewise,  no  measured  values  of  critical  point  data  are  available  for 

11  12  13 

any  metal  except  Hg  so  the  best  recent  estimates  were  used.  ’  ’ 

The  scaling  constants  are  summarized  in  Table  I. 


Table  I.  Scaling  Constants 


Metal 

RTC 

(RTr 

°K 

°K 

kJ/gm 

m/  sec 

Al 

933 

8650 

2.67 

1633 

Cu 

1356 

8500 

1.09 

1043 

Ag 

1234 

7460 

0.575 

758 

Au 

1336 

9500 

0.401 

630 

Hg 

234 

1733 

O.0718 

268 

Fb 

601 

5400 

0.217 

465 

The  theoretical  curve  is  computed  from  a  modified  van  der  Waals 
model,  in  dimensionless  form,  with  critical  compressibility  3/8?  and 
with  Cy(liq)  =  (5/2)R  and  Cv (vap)  =  (3/2)R  .  This  choice  of 
specific  heats,  made  possible  by  relinquishing  Maxwell's  equal  area 
rule,  is  a  better  approximation  to  the  known  specific  heat  data  than 
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assumption  of  identical  specific  heats  for  liquid  and  vapor.  In  its 
scaled  form  c  is  then  a  function  of  dimensionless  variables  only. 

The  scale  of  reduced  specific  heat  content,  q,  for  the  theoret¬ 
ical  curve  is  that  computed  from  the  van  der  Waals  function,  with  an 
added  0.1  unit  to  adjust  the  qM  value  so  that  it  represents,  reason¬ 
ably  well,  the  corresponding  values  for  the  six  metals  which  cluster 
fairly  closely  to  the  nominal,  empirical  value  of  0.60 

As  may  be  seen  in  Figure  3  the  data  and  the  theoretical  curve  • 
agree  quite  well  on  the  point  of  inception  of  vaporization  wave 
speed,  and  on  the  early  rise  of  the  curve.  The  data  points  for  the 
noble  metals  fall  below  the  theoretical  curve  and  go  out  beyond  the 

critical  point  limit  toward  what  appears  to  be  a  horizontal  asymptote. 

3 

Handbook  wave  speeds  for  Cu  actually  agree  better  with  the  plotted 
data  than  does  the  van  der  Waals  curve.  Some  recent  unpublished  data 
on  Cu  taken  in  our  laboratory,  reduced  by  a  finer-grained  method  of 
reduction,  adhere  more  closely  to  the  van  der  Waals  curve  up  to  an 
ordinate  of  0.15  before  bending  away.  If  the  correction  for  thermal 
expansion  could  be  applied,  the  wave  speed  values  would  be  increased 
by  at  least  20  -  70$>  depending  on  q,  the  larger  increases  applying 
near  critical.  Such  a  correction  would  somewhat  improve  the  agree¬ 
ment  for  Ag,  Cu,  Au  and  Al;  but  would  move  Fb  and  Hg  points  further 
away  from  the  theoretical  curve. 

Some  especial  problems  were  encountered  with  the  Hg  wires  which 
had  to  be  frozen  in  an  acetone  dry- ice  mixture,  kept  in  a  cold  box 
and  transferred  rapidly  to  the  test  cell  for  explosion.  Refinements 
in  this  process  might  result  in  better  reproducibility  of  the  mercury 

U.cl>  O  CL  . 
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8.  DEVIATIONS  FROM  THE  MODEL 

If  we  confine  our  attention  to  phenomena  below  critical*  two 
deviations  from  the  assumed  model  appear  to  be  of  primary  importance. 

The  theoretical  wave  speed  applies  to  the  head  of  an  expansion 
wave  traveling  into  the  molten  fluid.  Of  necessity  the  front  across 
which  the  conductivity  drops  to  zero  must  come  somewhere  in  the 
region  of  reduced  density  behind  the  head  of  the  wave.  It  may  occur 
across  a  narrow  segment  of  the  expansion  fan  and  thus  be  representable 
by  a  sharp  front  possibly  resembling  a  ’’conduct ivity  shock  wave,”  but 
its  speed  of  propagation  will  be  less  than  that  of  the  sound  wave 
which  first  travels  into  the  interior.  Because  we  know  little  about 
the  change  of  conductivity  in  the  expansion  wave  we  have  no  way  of 
estimating  the  decrement  in  wave  speeds  to  be  expected.  The  noble 
metals  and  A1  fall  below  the  theoretical  curve  in  the  expected  way 
but  Fb  and  Hg  do  not. 

It 

The  x-ray  pictures  of  Fansler  and  Shear  demonstrate  that  for 
wires  2-5  times  larger  and  for  somewhat  slower  rates  of  energy  ad¬ 
dition  the  expanding  cylinder  of  molten  fluid  is  not  uniform  but  dis¬ 
plays  the  transverse  variations  in  density  called  striations.  How 
these  striations  are  related  to  the  expansion  wave  is  not  known  at 
present  although  some  speculative  hypotheses  can  be  advanced.  Their 
presence  in  these  cases  appears  to  rule  out  the  hypothesis  that  the 
vaporization  wave  always  proceeds  uniformly  through  a  cylinder  of 
uniform  density. 

The  effect  of  unvaporized  portions  of  the  wire  remaining  in 
striations  appears  to  move  all  data  points  toward  higher  values  of  n 
by  variable  amounts.  There  is  presently  no  method  of  estimating 
this  effect.  Experiments  now  underway  are  directed  toward  determin¬ 
ing  the  relationship  between  number  of  striations  and  rate  of  energy 
addition  during  vaporization  wave  passage.  The  present  x-ray 
techniques  are  unable  to  resolve  any  striations  in  the  smaller  wires 
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and  faster  explosions  typical  of  much  of  the  data  plotted  in  Figure  3 
however 5  their  presence  in  the  larger,  slower  cases  is  a  cautionary 
warning  that  density  variations  somehow  play  an  important  part  in  the 
heating  and  expansion  of  superheated  metals. 
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Figure  1.  Sketch  of  two-phase  region  and  isotherms. 
P  -  pressure,  V  -  specific  volume. 
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Figure  2.  Streak  photograph  of  expanding  copper  wire  correlated  in 
time  with  current  and  voltage, 
d  =  .0254  cm,  V0  =  3kV,  C  =  32  p,F. 
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Figure  3 •  Wave  speeds  for  several  metals. 
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satisfactory  but  higher  values  deviate  considerably  from  theory.  Possible  causes  of 
the  deviations  include  the  crudity  of  the  fluid  dynamic  model,  neglect  of  thermal  ex¬ 
pansion,  lack  of  information  about  the  relationship  between  density  and  electrical 
conductivity  and  the  approximation  imposed  by  the  van  der  Waals  equation. 


4  A  RIPLAClt  DO  FORM  MTS.  t  JAN  M,  WHICH  IB 

1  A  J  A  OBSOLETE  FOR  ARMY  UBS, 


UNCLASSIFIED 


curlty  Classification 


key  wo  ft  os 


Vaporization  Waves 
Fluid  Dynamics 
High  Temperature  Metals 
Unsteady  Flows 
Simple  Waves 
Critical  Phenomena 


UNCLASSIFIED 

Security  Classification 


